Long term surface soil moisture (SSM) data with stable and consistent quality are critical for global environment and climate change monitoring. L band radiometers onboard the recently lunched Soil Moisture Active Passive (SMAP) Mission can provide the state-of-the-art accuracy SSM, but the short temporal coverage of the data records has limited its applications in long-term studies. While Advanced Microwave Scanning Radiometer for EOS (AMSR-E) and AMSR2 series provide long term observational records of multi-frequency radiometers (C, X, and K bands). This dataset contains 20 years (2002-2022) global spatio-temporal consistent surface soil moisture. The resolution is 36 km at daily scale, the projection is EASE-Grid2, and the data unit is m3 / m3. This dataset adopts the soil moisture neural network retrieval algorithm developed by Yao et al. (2017). This study transfers the merits of SMAP to AMSR-E/2 through using an Artificial Neural Network (ANN) in which SMAP standard SSM products serve as training targets with AMSR-E/2 brightness temperature (TB) as input. Finally, long term soil moisture data are output. This dataset can reproduce the spatial and temporal distribution of SMAP soil moisture, with comparable accuracy as SMAP soil moisture product. This dataset also compares well with in situ SSM observations at 14 dense validation networks globally, with accuracy of 5% volumetric water content, and outperforms AMSR-E/2 standard SSM products from JAXA and LPRM. This global observation-driven dataset spans nearly two decades at present, and is extendable though the ongoing AMSR2 and upcoming AMSR3 missions for long-term studies of climate extremes, trends, and decadal variability.
YAO Panpan, LU Hui
This dataset contains daily 0.01°×0.01° land surface soil moisture products in the Qinghai-Tibet Plateau in 2005, 2010, 2015, 2017, and 2018. The dataset was produced by utilizing the multivariate statistical regression model to downscale the “SMAP Time-Expanded 0.25°×0.25° Land Surface Soil Moisture Dataset in the Qinghai-Tibet Plateau (SMsmapTE, V1)”. The auxiliary datasets participating in the multivariate statistical regression include GLASS Albedo/LAI/FVC, 1km all-weather surface temperature data in western China by Ji Zhou, and Lat/Lon information.
CHAI Linna, ZHU Zhongli, LIU Shaomin
Based on 11 well-acknowledged global-scale microwave remote sensing-based surface soil moisture products, and with 9 main quality impact factors of microwave-based soil moisture retrieval incorporated, we developed the Remote Sensing-based global Surface Soil Moisture dataset (RSSSM, 2003~2020) through a complicated neural network approach. The spatial resolution of RSSSM is 0.1°, while the temporal resolution is approximately 10 days. The original dataset covered 2003~2018, but now it has been updated to 2020. RSSSM dataset is outstanding in terms of temporal continuity, and has full spatial coverage except for snow, ice and water bodies. The comparison against the global-scale in-situ soil moisture measurements indicates that RSSSM has a higher spatial and temporal accuracy than most of the frequently-used global/regional long-term surface soil moisture datasets. In addition, although RSSSM is remote sensing based, without the incorporation of any precipitation data or records, its interannual variation generally conforms with that of precipitation (e.g., the GPM IMERG precipitation data) and Standardized Precipitation Evapotranspiration Index (SPEI). Moreover, RSSSM can also reflect the impact of human activities, e.g., urbanization, cropland irrigation and afforestation on soil moisture changes to some degree. The data is in ‘Tiff’ format, and the size after compression is 2.48 GB. The relevant data describing paper has been published in the Journal ‘Earth System Science Data’ in 2021.
CHEN Yongzhe, FENG Xiaoming, FU Bojie
The long-term evolution of lakes on the Tibetan Plateau (TP) could be observed from Landsat series of satellite data since the 1970s. However, the seasonal cycles of lakes on the TP have received little attention due to high cloud contamination of the commonly-used optical images. In this study, for the first time, the seasonal cycle of lakes on the TP were detected using Sentinel-1 Synthetic Aperture Radar (SAR) data with a high repeat cycle. A total of approximately 6000 Level-1 scenes were obtained that covered all large lakes (> 50 km2) in the study area. The images were extracted from stripmap (SM) and interferometric wide swath (IW) modes that had a pixel spacing of 40 m in the range and azimuth directions. The lake boundaries extracted from Sentinel-1 data using the algorithm developed in this study were in good agreement with in-situ measurements of lake shoreline, lake outlines delineated from the corresponding Landsat images in 2015 and lake levels for Qinghai Lake. Upon analysis, it was found that the seasonal cycles of lakes exhibited drastically different patterns across the TP. For example, large size lakes (> 100 km2) reached their peaks in August−September while lakes with areas of 50−100 km2 reached their peaks in early June−July. The peaks of seasonal cycles for endorheic lakes were more pronounced than those for exorheic lakes with flat peaks, and glacier-fed lakes with additional supplies of water exhibited delayed peaks in their seasonal cycles relative to those of non-glacier-fed lakes. Large-scale atmospheric circulation systems, such as the westerlies, Indian summer monsoon, transition in between, and East Asian summer monsoon, were also found to affect the seasonal cycles of lakes. The results of this study suggest that Sentinel-1 SAR data are a powerful tool that can be used to fill gaps in intra-annual lake observations.
ZHANG Yu, ZHANG Guoqing
Agricultural irrigation consumes a large amount of available freshwater resources and is the most immediate human disturbance to the natural water cycle process, with accelerated regional water cycles accompanied by cooling effects. Therefore, estimating irrigation water use (IWU) is important for exploring the impact of human activities on the natural water cycle, quantifying water resources budget, and optimizing agricultural water management. However, the current irrigation data are mainly based on the survey statistics, which is scattered and lacks uniformity, and cannot meet the demand for estimating the spatial and temporal changes of IWU. The Global Irrigation Water Use Estimation Dataset (2011-2018) is calculated by the satellite soil moisture, precipitation, vegetation index, and meteorological data (such as incoming radiation and temperature) based on the principle of soil water balance. The framework of IWU estimation in this study coupled the remotely sensed evapotranspiration process module and the data-model fusion algorithm based on differential evolution. The IWU estimates provided from this dataset have small bias at different spatial scales (e.g., regional, state/province and national) compared to traditional discrete survey statistics, such as at Chinese provinces for 2015 (bias = −3.10 km^3), at U.S. states for 2013 (bias = −0.42 km^3), and at various FAO countries (bias = −10.84 km^3). Also, the ensemble IWU estimates show lower uncertainty compared to the results derived from individual precipitation and soil moisture satellite products. The dataset is unified using a global geographic latitude and longitude grid, with associated metadata stored in corresponding NetCDF file. The spatial resolution is about 25 km, the time resolution is monthly, and the time span is 2011-2018. This dataset will help to quantitatively assess the spatial and temporal patterns of agricultural irrigation water use during the historical period and support scientific agricultural water management.
ZHANG Kun, LI Xin, ZHENG Donghai, ZHANG Ling, ZHU Gaofeng
Surface soil moisture (SSM) is a crucial parameter for understanding the hydrological process of our earth surface. Passive microwave (PM) technique has long been the primary choice for estimating SSM at satellite remote sensing scales, while on the other hand, the coarse resolution (usually >~10 km) of PM observations hampers its applications at finer scales. Although quantitative studies have been proposed for downscaling satellite PM-based SSM, very few products have been available to public that meet the qualification of 1-km resolution and daily revisit cycles under all-weather conditions. In this study, therefore, we have developed one such SSM product in China with all these characteristics. The product was generated through downscaling of AMSR-E and AMSR-2 based SSM at 36-km, covering all on-orbit time of the two radiometers during 2003-2019. MODIS optical reflectance data and daily thermal infrared land surface temperature (LST) that have been gap-filled for cloudy conditions were the primary data inputs of the downscaling model, in order to achieve the “all-weather” quality for the SSM downscaling outcome. Daily images from this developed SSM product have achieved quasi-complete coverage over the country during April-September. For other months, the national coverage percentage of the developed product is also greatly improved against the original daily PM observations. We evaluated the product against in situ soil moisture measurements from over 2000 professional meteorological and soil moisture observation stations, and found the accuracy of the product is stable for all weathers from clear sky to cloudy conditions, with station averages of the unbiased RMSE ranging from 0.053 vol to 0.056 vol. Moreover, the evaluation results also show that the developed product distinctly outperforms the widely known SMAP-Sentinel (Active-Passive microwave) combined SSM product at 1-km resolution. This indicates potential important benefits that can be brought by our developed product, on improvement of futural investigations related to hydrological processes, agricultural industry, water resource and environment management.
SONG Peilin, ZHANG Yongqiang
1. Glacial lake data sets (1960s−2020) This data set contains glacial lake data for the 1960s, 2016, 2017, 2018, 2019, and 2020, mapped from Korona KH-4, Sentinel-2, and Sentinel-1 imagery. 2. Potential Outburst Flood Hazard level of Bhutanese glacial lakes This data contains the Potential Outburst Flood Hazard level of Bhutanese glacial lakes with an area greater than 0.05 km2 (n=278). The value for each hazard assessment criteria is also provided in the data attributes.
RINZIN Sonam, ZHANG Guoqing
The Antarctic ice sheet elevation data were generated from radar altimeter data (Envisat RA-2) and lidar data (ICESat/GLAS). To improve the accuracy of the ICESat/GLAS data, five different quality control indicators were used to process the GLAS data, filtering out 8.36% unqualified data. These five quality control indicators were used to eliminate satellite location error, atmospheric forward scattering, saturation and cloud effects. At the same time, dry and wet tropospheric, correction, solid tide and extreme tide corrections were performed on the Envisat RA-2 data. For the two different elevation data, an elevation relative correction method based on the geometric intersection of Envisat RA-2 and GLAS data spot footprints was proposed, which was used to analyze the point pairs of GLAS footprints and Envisat RA-2 data center points, establish the correlation between the height difference of these intersection points (GLAS-RA-2) and the roughness of the terrain relief, and perform the relative correction of the Envisat RA-2 data to the point pairs with stable correlation. By analyzing the altimetry density in different areas of the Antarctic ice sheet, the final DEM resolution was determined to be 1000 meters. Considering the differences between the Prydz Bay and the inland regions of the Antarctic, the Antarctic ice sheet was divided into 16 sections. The best interpolation model and parameters were determined by semivariogram analysis, and the Antarctic ice sheet elevation data with a resolution of 1000 meters were generated by the Kriging interpolation method. The new Antarctic DEM was verified by two kinds of airborne lidar data and GPS data measured by multiple Antarctic expeditions of China. The results showed that the differences between the new DEM and the measured data ranged from 3.21 to 27.84 meters, and the error distribution was closely related to the slope.
HUANG Huabin
This biophysical permafrost zonation map was produced using a rule-based GIS model that integrated a new permafrost extent, climate conditions, vegetation structure, soil and topographic conditions, as well as a yedoma map. Different from the previous maps, permafrost in this map is classified into five types: climate-driven, climate-driven/ecosystem-modified, climate-driven/ecosystem protected, ecosystem-driven, and ecosystem-protected. Excluding glaciers and lakes, the areas of these five types in the Northern Hemisphere are 3.66×106 km2, 8.06×106 km2, 0.62×106 km2, 5.79×106 km2, and 1.63×106 km2, respectively. 81% of the permafrost regions in the Northern Hemisphere are modified, driven, or protected by ecosystems, indicating the dominant role of ecosystems in permafrost stability in the Northern Hemisphere. Permafrost driven solely by climate occupies 19% of permafrost regions, mainly in High Arctic and high mountains areas, such as the Qinghai-Tibet Plateau.
RAN Youhua, M. Torre Jorgenson, LI Xin, JIN Huijun, Wu Tonghua, Li Ren, CHENG Guodong
This dataset contains daily 0.05°×0.05° land surface soil moisture products in Qilian Mountain Area in 2019. The dataset was produced by utilizing the optimized wavelet-coupled-RF downscaling model (RF-OWCM) to downscale the “AMSR-E and AMSR2 TB-based SMAP Time-Expanded Daily 0.25°×0.25° Land Surface Soil Moisture Dataset in Qilian Mountain Area (SMsmapTE, V1)”. The auxiliary datasets participating in the downscaling model include GLASS Albedo/LAI/FVC, Thermal and Reanalysis Integrating Medium-resolution Spatial-seamless LST – Tibetan Plateau (TRIMS LST-TP) by Ji Zhou and Lat/Lon information.
CHAI Linna, ZHU Zhongli, LIU Shaomin
The Qinghai-Tibetan Plateau (QTP), the largest high-altitude and low-latitude permafrost zone in the world, has experienced rapid permafrost degradation in recent decades, and one of the most remarkable resulting characteristics is the formation of thermokarst lakes. Such lakes have attracted significant attention because of their ability to regulate carbon cycle, water, and energy fluxes. However, the distribution of thermokarst lakes in this area remains largely unknown, hindering our understanding of the response of permafrost and its carbon feedback to climate change.Based on more than 200 sentinel-2A images and combined with ArcGIS, NDWI and Google Earth Engine platform, this data set extracted the boundary of thermokarst lakes in permafrost regions of the Qinghai-Tibet Plateau through GEE automatic extraction and manual visual interpretation.In 2018, there were 121,758 thermokarst lakes in the permafrost area of the Qinghai-Tibet Plateau, covering an area of 0.0004-0.5km², with a total area of 1,730.34km² respectively.The cataloging data set of Thermokarst Lakes provides basic data for water resources evaluation, permafrost degradation evaluation and thermal karst study on the Qinghai-Tibet Plateau.
CHEN Xu, MU Cuicui, JIA Lin, LI Zhilong, FAN Chengyan, MU Mei, PENG Xiaoqing, WU Xiaodong
This dataset contains land surface soil moisture products with SMAP time-expanded daily 0.25°×0.25°in Qinghai-Tibet Plateau Area. The dataset was produced based on the Random Forest method by utilizing passive microwave brightness temperature along with some auxiliary datasets. The temporal resolution of the product in 1980,1985,1990,1995 and 2000 is monthly, by using SMMR, SSM/I, and SSMIS brightness temperature from 19 GHz V/H and 37 GHz V channels. The temporal resolution of the product between June 20, 2002 and Dec 30, 2018 is daily, by utilizing AMSR-E and AMSR2 brightness temperature from 6.925 GHz V/H, 10.65 GHz V/H, and 36.5 GHz V channels. The auxiliary datasets participating in the Random Forest training include the IGBP land cover type, GTOPO30 DEM, and Lat/Lon information.
CHAI Linna, ZHU Zhongli, LIU Shaomin
1、Based on field eddy correlation (EC) measurement data, using the standard data processing method for EC data, including despiking, coordinate rotation, air density corrections, outlier rejection, and friction velocity threshold (u*) corrections, gap filled, and NEE partition. The dataset collects carbon flux data and microclimate measurement data from 2003 to 2016 in three typical alpine grassland ecosystems on the Qinghai-Tibet Plateau, including Damxung alpine meadow, Haibei alpine meadow ，Naqu alpine meadow，Zoige alpine grassland，Qilian mountion grassland . The time resolution of data is high (30 min), and the interpolation of data is complete throughout the year. This dataset can be applied to carbon flux assessment, comparison and prediction in these alpine meadows, attribution of climate factors affecting carbon flux, validation of model simulation results, etc. 2、Based on the MCDGF43 dataset, we produce the visible and near-infared albedo of Tibetan Plateau, using the standard data processing of hdf to tif , including the moasic, resample and masked by Tibetan Plateau's boundary. The time resolution of dataset is 8 days and the spatial resolution is 500 meters, which span the period of 2003-2016.
ZHANG Yangjian, SU Peixi, YANG Yan
The proportion data set of daily cloudless MODIS snow cover area in babaohe river basin (2008.1.1-2014.6.1) was obtained after cloud removal processing using a cloud removal algorithm based on cubic spline function interpolation on the basis of daily cloudless MODIS snow cover product-mod10a1 (tang zhiguang, 2013). This data set adopts the projection method of UTM (horizontal axis isometric cutting cylinder), with a spatial resolution of 500m, and provides Daily Snow Albedo daily-sad results for the babao river basin.The data set is a daily file from January 1, 2008 to June 1, 2014.Each file is the snow albedo result of the day, with a value of 0-100 (%), is the ENVI standard file, and the naming rule is: mod10a1.ayyyyddd_h25v05_snow_sad_grid_2d_reproj_babaohe_nocloud.img, where YYYY represents the year, DDD stands for Julian day (001-365/366).The file can be opened directly with ENVI or ARCMAP software. The original MODIS snow cover data products processed by declouding are derived from MOD10A1 products processed by the us national snow and ice data center (NSIDC). This data set is in HDF format and USES sinusoidal projection. The attributes of the cloud-free MODIS albedo data set (2008.1.1-2014.1.1) in babaohe river basin are composed of the spatial and temporal resolution, projection information and data format of the dataset.
WANG Jian, PAN Haizhu
This dataset contains SMAP time-expanded monthly 0.25°×0.25° land surface soil moisture products in Qilian Mountain Area of 1980, 1985, 1990, 1995 and 2000. The dataset was produced based on Random Forest method by utilizing SMMR, SSM/I and SSMIS brightness temperature from 19 GHz V/H and 37 GHz V channels along with some auxiliary datasets. The auxiliary datasets participating in the RF training include IGBP land cover type, GTOPO30 DEM and Lat/Lon information.
CHAI Linna, ZHU Zhongli, LIU Shaomin
The dataset of ground truth measurement synchronizing with Envisat ASAR was obtained in the arid region hydrological experimental area on Sep. 19, 2007 during the pre-observation period. One scene of Envisat ASAR image was captured on Sep. 19. The data were in AP mode and VV/VH polarization combinations, and the overpass time was approximately at 11:29 BJT. Those provide reliable ground data for remote sensing retrieval and validation of soil moisture from Envisat ASAR image. Observation items included: (1) soil moisture measured by the cutting ring method in Linze reed land, Zhangye farmland, Zhangye gobi, Linze maize land, Linze alfalfa land, Zhangye weather station, and Linze wetland. (2) GPS measured by GARMIN GPS 76 (3) vegetation measurements including the vegetation height, the green weight, the dry weight, the sampling method, and descriptions on the land type, uniformity and dry and wet conditions (4) atmospheric parameters at Daman Water Management office measured by CE318 (produced by CIMEL in France). The total optical depth, aerosol optical depth, Rayleigh scattering coefficient, column water vapor in 936 nm, particle size spectrum and phase function were then retrieved from these observations. The optical depth in 1020nm, 936nm, 870nm, 670nm and 440nm were all acquired by CE318. Those data include the raw data in .k7 and can be opened by ASTPWin. ReadMetext files (.txt) is attached for detail. Processed data (after retrieval of the raw data) archived as Excel files are on optical depth, rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number. (5) roughness measured by the roughness plate together with the digital camera. The coordinates of the sample would be got with the help of ArcView; and after geometric correction, surface height standard deviation (cm) and correlation length (cm) could be acquired based on the formula listed on pages 234-236, Microwave Remote Sensing (Vol. II). The roughness data were initialized by the sample name, which was followed by the serial number, the name of the file, standard deviation and correlation length. Each text files (.txt) file is matched with one sample photo and standard deviation and correlation length represent the roughness. In addition, the length of 101 radius is also included for further checking.
CHE Tao, LI Xin, BAI Yunjie, DING Songchuang, GAO Song, HAN Xujun, HAO Xiaohua, LI Hongyi, LI Zhe, LIANG Ji, PAN Xiaoduo, QIN Chun, RAN Youhua, WANG Xufeng, WU Yueru, YAN Qiaodi, ZHANG Lingmei, FANG Li, LI Hua, Liu Qiang, Wen Jianguang, MA Hongwei, YAN Yeqing, YUAN Xiaolong
The dataset of ground truth measurement synchronizing with ALOS PALSAR was obtained in the Linze grassland foci experimental area on Jun. 10, 2008. The data were in FBS mode and HH/HV polarization combinations, and the overpass time was approximately at 23:39 BJT. Observations were carried out in plots A, B, C, D and E, which were divided into 6×6 subsites, with each one spanning a 120×120 m2 plot. Soil gravimetric moisture, volumetric moisture, and soil bulk density after drying by the cutting ring and the mean soil temperature from 0-5cm by the probe thermometer were measured in A, B and C; the soil temperature, soil moisture, the loss tangent, soil conductivity, the real part and the imaginary part of soil complex permittivity by the POGO soil sensor, and the mean soil temperature from 0-5cm by the probe thermometer in D and E. Data were archived in Excel file. See WATER: Dataset of setting of the sampling plots and stripes in the foci experimental area of Linze station for more information.
BAI Yanfen, CAO Yongpan, GE Chunmei, HU Xiaoli, WANG Shuguo, Wang Weizhen, WU Yueru, ZHU Shijie, FENG Lei
The dataset of ground truth measurements synchronizing with Envisat ASAR was obtained in C1, W2 and B2 of the Biandukou foci experimental area on Mar. 14, 2008, from 23:30 on 14 to 1:00 on 15, to be specific. The ASAR data were in AP mode and VV/VH polarization combinations, and the overpass time was approximately at 23:21 BJT. The wheat stubble land, the deep plowed land and the rape stubble land were chosen for measurements. (1) The surface radiative temperature and the physical temperature were measured by the handheld infrared thermometer. Besides, the land cover type was also recorded. The data can be opened by Microsoft Office. (2) The gravimetric soil moisture (samples from 0-1cm, 1-3cm, 3-5cm, 5-10cm and 10-20cm) was measured by the microwave drying method. (3) The frost depth by the chopstick and the ruler. The soil was considered frozen when it was hard and with ice crystal. The data can be opened by Microsoft Office. Four data files were included, ASAR data, C1, W2 and B2 data.
CHANG Sheng, Fang Qian, QU Ying, LIANG Xingtao, LIU Zhigang, PAN Jinmei, PENG Danqing, REN Huazhong, ZHANG Yongpan, ZHANG Zhiyu, ZHAO Shaojie, Zhao Tianjie, ZHENG Yue, Zhou Ji, LIU Chenzhou, YIN Xiaojun, ZHANG Zhiyu
The dataset of ground truth measurements synchronizing with Envisat ASAR was obtained in No. 1 and 2 quadrates of the Biandukou foci experimental area on Oct. 18, 2007, during the pre-observation period. The ASAR data were in AP mode and VV/VH polarization combinations, and the overpass time was approximately at 11:17 BJT. Both the quadrates were divided into 3×3 subsites, with each one spanning a 30×30 m2 plot. 25 sampling points were chosen, including centers and corners. Simultaneous with the satellite overpass, numerous ground data were collected: the soil temperature , volumetric soil moisture (cm^3/cm^3), soil salinity (s/m), soil conductivity (s/m) by the Hydra probe, the surface radiative temperature by the handheld infrared thermometer, gravimetric soil moisture, volumetric soil moisture, and soil bulk density by drying soil samples from the cutting ring (100cm^3). Meanwhile, vegetation parameters as height, coverage and water content were also observed. Those provide reliable ground data for the development and validation of soil moisture, soil freeze/thaw algorithms and the forward model from active remote sensing approaches.
BAI Yunjie, CAO Yongpan, WANG Jian, Wang Weizhen, WANG Xufeng, JIN Rui, Qu Yonghua, ZHOU Hongmin
The dataset of ground truth measurements synchronizing with Envisat ASAR was obtained in the Linze grassland foci experimental area on Jul. 11, 2008. The data were in AP mode and VV/VH polarization combinations, and the overpass time was approximately at 11:26 BJT. Observations were carried out in the reed plot A, the saline plots B and C, the alfalfa plot D and the barley plot E, which were divided into 6×6 subsites, with each one spanning a 120×120 m2 plot. Soil gravimetric moisture, volumetric moisture, and soil bulk density after drying by using the cutting ring, the mean soil temperature from 0-5cm by the probe thermometer, and the canopy temperature and the land surface temperature by the hand-held infrared thermometer were measured in A, B and C; the soil temperature, soil moisture, the loss tangent, soil conductivity, the real part and the imaginary part of soil complex permittivity by the POGO soil sensor, the mean soil temperature from 0-5cm by the probe thermometer, the canopy temperature and the land surface temperature by the hand-held infrared thermometer in D and E. Data were archived in Excel file. See WATER: Dataset of setting of the sampling plots and stripes in the foci experimental area of Linze station for more information.
CAO Yongpan, CHAO Zhenhua, GE Chunmei, HU Xiaoli, HUANG Chunlin, LIU Chao, WU Yueru, SHEN Xinyi
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS
0931-4967287
poles@itpcas.ac.cn
Links
National Tibetan Plateau Data CenterFollow Us
A Big Earth Data Platform for Three Poles © 2018-2020 No.05000491 | All Rights Reserved
｜ 
No.11010502040845
Tech Support: westdc.cn
